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ABSTRACT 

Using NASA IRTF SpeX data from 0.8 to 4.5 /zm, we determine self-consistently the stellar prop- 
erties and excess emission above the photosphere for a sample of classical T Tauri stars (CTTS) in 
the Taurus molecular cloud with varying degrees of accretion. This process uses a combination of 
techniques from the recent literature as well as observations of weak-line T Tauri stars (WTTS) to 
account for the differences in surface gravity and chromospheric activity between the TTS and dwarfs, 
which are typically used as photospheric templates for CTTS. Our improved veiling and extinction 
estimates for our targets allow us to extract flux-calibrated spectra of the excess in the near-infrared. 
We find that we are able to produce an acceptable parametric fit to the near-infrared excesses using 
a combination of up to three blackbodies. In half of our sample, two blackbodies at temperatures 
of 8000 K and 1600 K suffice. These temperatures and the corresponding solid angles are consistent 
with emission from the accretion shock on the stellar surface and the inner dust sublimation rim of 
the disk, respectively. In contrast, the other half requires three blackbodies at 8000, 1800, and 800 K, 
to describe the excess. We interpret the combined two cooler blackbodies as the dust sublimation wall 
with either a contribution from the disk surface beyond the wall or curvature of the wall itself, neither 
of which should have single-temperature blackbody emission. In these fits, we find no evidence of a 
contribution from optically thick gas inside the inner dust rim. 

Subject headings: open cluster and associations: individual (Taurus) — stars: pre-main sequence — 
infrared: stars 



1. INTRODUCTION 

It is now accepted, from studies of UV and optical spec- 
troscopy and near-infrared (NIR) interferometry, that 
the excess emission shortward of 5 /xm in spectral energy 
distributions (SEDs) of classical T Tauri stars (CTTS) 
arises mainly in two physical components. The first, re- 
sponsible for the near-infrared (NIR) excess, is the sharp 
inner edge, or 'wall', of the dust disk at which the disk 
temperature becom es high enough to sublimate the dust 
(|Natta et al.ll200ll ). Inside the dust-sublimation radius 
only gas remains, spiraling in until it is accreted onto 
the star through magnetospheric accretion columns. The 
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second emission component, responsible for the UV and 
optical excess, is the accretion shock formed when mate- 
rial free-falling along stellar magnetic field lines merges 
with the stellar photosphere through a shock at the sur- 
face, where most of the a ccretion luminosity is emitted 
(|Calvet fc GullbrindH998l) . 

One can constrain the accretion rate and the struc- 
ture of the inner disk wall from the shape, absolute 
flux, and line emission seen in the excess. However, 
extracting the excess from the observed spectrum re- 
quires subtracting the underlying stellar photosphere 
and correcting for extinction along the line of sight. 
This is typically accomplished by identifying the pho- 
tospheric template that best matches the intrinsic stel- 
lar spectrum of the program star, assuming that the 
weaker absorption features in the program star are due 
to the excess continuum emission filling in, or 'veil- 
ing', these features. By comparing the difference in 
feature depth or equivalent width between the program 
star and template it is possible to determine the frac- 
tion of th e observed flux contribu t ed by the stellar pho- 
tosphcrc (Basri & Batalha 1990; Harti gan et al.l 119951 : 
IMuzerolle et al.ll2003t lEspaillat et al.ll2010D . The slopes 
of the veiling-corrected spectrum over a range of wave- 
lengths can be compared with that of the intrinsic tem- 
plate to derive the extinction, Ay. Subtraction of the 
intrinsic template, offset from the observed spectrum by 
the appropriate degree of veiling, from the extinction- 
corrected observed spect r um yields the veiling spectrum 
(iGullbring et all Il998at iSeperuelo Duarte et al.l 120081: 
iFischer et al.l I2011IL If the derived excess spectrum is 
absolutely flux calibrated, e.g. with simultaneous pho- 
tometry, then one also can obtain sizes of its emitting 
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regions. 

The NIR excess is of particular interest as it mea- 
sures directly the inner disk emission, thus providing in- 
sight into the state of the inner disk and its effect on 
the outer disk. For example, the maximum grain size 
and composition in the wall dust population affect its 
shape and hence the emi tting area (jlsella k Nattal [20051 : 
iTannirkulam et aT]|2007f l. There is also evidence for an 
optically th ick gas component inside the dust sublima- 
tion radius (ITannirkulam et all 120081 : lEisner et al.ll2007t 
iFischer et al. 2011, and references therein). The geom- 
etry of the wall relative to the disk behind it deter- 
mines how much stellar emission is incident on the outer 
disk; if the wall 'shadows' the disk it would prevent the 
middle few AU from being heated effectively, produc- 
ing less flaring (iNatta et alJl200l1lDullemond et al.ll200ll 
iMeeus et al1l200lHDullemond k Don iinik 2QM) 

In practice it can be difficult to successfully apply the 
standard approach for excess extraction. Careful selec- 
tion is required to avoid lines that show differential veil- 
ing due to chromospheric em ission (|Finkenzeller k Basril 
119871 Batalha k Basri 1993). It is also essential to se- 
lect the correct template star; the best choice is a weak- 
line T Tauri stars (WTTS), which does not accrete and 
has no NIR excess. These stars are quantitatively dis- 
tinguished fro m the CTTS by the st rength of their Ha 
emission lines ( White k Ghezl 1 20 01). Of possible tem- 
plate stars, WTTS have the closest physical properties 
to the CTTS, including active chromospheres, compara- 
ble metallicities, and similar surface gravities. The lat- 
ter is an essential consideration, as TTS have surface 
gravities between those of dwarfs and giants. Studies 
comparing these templates indicate that dwarf standard 
stars provide an acceptable match to TTS photospheres 
at optical wavelengths (Basri k Batalha 1990). How- 
ever, many of the strong NIR absorption features are 
gravity-sensitive, and the effect of surface gravity on the 
molecular bands near 1.65 fj,m produces a distinctive tri- 
angular shape in the continuum emission that must be 
accounte d for when determining the excess over that re- 
gion, (see [Gr^XSorbaIii][200| Fig. 9.6) 

Here we use an infrared spectrograph, SpeX, to study 
the NIR excess in ten accreting TTS. SpeX is the ideal 
instrument for this project, as it obtains continuous 0.8 
to 2.5 fim coverage. In this paper, we develop meth- 
ods of veiling determination that account for the inter- 
mediate surface gravities of TTS and identify a set of 
lines that appear to have minimal effects from chromo- 
spheric activity. Using these veiling values and WTTS 
template stars, we extract the excess above the photo- 
sphere in each of our targets. We then perform a model- 
independent black body analysis to place constraints on 
the temperatures and emitting areas of the excess com- 
ponent, testing in the process the possibility of optically 
thick emission from inside the dust sublimation radius. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection of CTTS and WTTS 

Our initial sample of eight stars was selected from 
CTTS in the Taurus-Auriga molecular cloud complex 
for whic h we have Svitzer Infrared Spectrograph (IRS) 
spectra ([Fur Ian et al.l 120 06) that 1) do not indicate the 
presence of gaps or other radial structure in their disks, 



2) show a range of apparent excess in the NIR, based 
on 2MASS photometry and their spectral types in the 
literature, and 3) show a range of mass accretion rates 
from the literature. Additional observations of two bright 
WTTS were included at our two most common spec- 
tral types (K7 and M2) for use as photospheric tem- 
plates. The CTTS are all single stars, within a de- 
tection limit of A?n#-=2 at separations greater than 20 
mas, while the WT TS, LkCa 3 and V 827 Tau, are sub- 
arcsecond binaries (jKraus et al.ll20lH ). For the purposes 
of studying the excess over the photosphere, using bi- 
naries is obviously un-ideal. However, based on their 
colors, the companions to LkCa 3 and V827 Tau are 
M6 and M5 and contribute less than half of the to- 
tal flux in the NIR ( 45% at K R for LkCa 3 and 37% 
at H for V827 Tau) (White k Ghezl |200H IKraus et al.l 
120111) . Since these types will p e ak around H or K 
band (jKenvon k Hartmannl 119951: IRavner et all 120091) , 
they should be responsible for less flux than this at i, 
z, and J bands. Nonetheless, in Sj3] we demonstrate sev- 
eral techniques and checks to minimize the impact of the 
WTTS' binarity on our analysis. The 2MASS magni- 
tudes and stellar parameters of our sample are given in 
Tables [Hand [21 respectively. For each target we obtained 
spectra and photometry. 

2.2. Spectroscopy 

Th e spectra were obtained with SpeX (jRavner et al.1 
120031 ) at the NASA Infrared Telescope Facility (IRTF) 
on 1, 2, and 3 December 2010. We observed our targets 
with the long- wavelength, cross-dispersed 2.1 /im. mode 
(LXD) with the 0'.'5xl5'.'0 slit (R=A/AA=2000) cover- 
ing 2.1 to 4.5/im. Following each LXD observation, we 
also observed the targets in the short-wavelength, cross- 
dispersed mode (SXD) with the 0'.'3xl5'.'0 slit (R=2000) 
from 0.8 to 2.5 microns. Our SXD integration times were 
selected to have S/N > 100 in the H band, while our 
LXD times were selected for S/N ^20 at the K band. 
The data were obtained with the slit rotated to the par- 
allactic angle and in an ABBA nod pattern. 

The spectra were reduced with the Spextool package 
([Cushing et al.1 [20041 ). We sky-subtracted each individ- 
ual exposure using the opposite nod positions, extracted 
them separately, scaled each spectrum in the set to the 
collective median value, and combined them using the 
robust median option. Spextool includes a package to 
correct for telluric absorption and perform relative flux 
calibration (jVacca et al.l 120031 ) using observations of A0 
stars. These stars provide a relatively featureless contin- 
uum in the infrared against which the telluric absorption 
features can be clearly identified and removed. Relative 
flux calibration is achieved by comparison of the telluric 
standard with a high- resolution model of Vega. This pro- 
cess first requires the removal of the intrinsic hydrogen 
series absorption lines in the underlying spectrum of the 
telluric standard, which can be accomplished using ei- 
ther the instrument profile or a convolution between the 
telluric standard and Vega model. The advantage of con- 
volution is that it produces a better correction than the 
instrument profile; however, it requires an individual hy- 
drogen line with sufficiently high S/N. For these data, 
we were able to perform the convolution with Vega using 
the Paschen S and Brackett 7 lines, respectively. 

Unfortunately, half of our observations in SXD had 
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airmass differences greater than 0.1 between the targets 
and the nearest telluric standard observation. All but 
two of the target observations were bracketed in airmass 
between two telluric observations. For the targets with a 
poor match in airmass to any individual A0 star observa- 
tion, we were able to produce an improved correction by 
taking an average of the two spectra of the program star, 
each corrected with one of the bracketing telluric obser- 
vations and weighting by the airmass difference with the 
program star. 

Finally, the telluric-corrected, combined spectrum for 
a given order was merged with the neighboring order to 
produce a continuous spectrum. We then excised the 
portions of each spectrum corresponding to the large tel- 
luric bands with less than 10% transmission, e.g. 1.35 < 
A < 1.45, as these regions were too noisy to use. The 
LXD spectra were scaled to the SXD spectra between 
2.29 and 2.35 microns and spliced together at the point 
of equal S/N. The final spectra are displayed in Figure 
ffl 

2.3. Photometry 

To flux calibrate our spectra, following each observa- 
tion we imaged the target with the guide camera (0.1185 
arcsec/pixel) at if in a 7-position dither pattern. We 
then observed a photometric standard fro m the extende d 
list of UKIRT faint standards (H awarden et al.M200lD . 
which was typically fainter than the targets. Conditions 
were photometric, and integration times were chosen to 
maintain the same number of counts per pixel between 
the target and standard to avoid differences in the level of 
non-linearity. The binary WTTS were not resolved, and 
the photometry listed for them in Table [1] is for both 
components. The unregistered images of the photomet- 
ric standard were median combined to create a master 
flat, as they typically had longer exposure times. After 
applying the flat field image to the data for the target 
and standard, we registered and median-combined the 
flat-fielded images for each target. Photometry was ex- 
tracted from the final image using IRAF's phot routine, 
with a gain of 14 e~ /DN, an instrumental zero-point 
magnitude of 20.57 at K, an aperture of 10 pixels, and a 
sky annulus from 30 to 35 pixels. For each night we then 
derived a photometric solution using the observations of 
our standards at all airmasses and applied it to the data. 
The resulting photometry is given in the last column of 
Table [TJ roughly half of the targets differ by at most 0.1 
mag between our measurement and their 2MASS magni- 
tudes, while the other half differ by up to 0.5 mag. Since 
T Tauri stars are highly variable, these differences are 
likely genuine. 

2.4. Ancillary Spectroscopy 

To test our infrared spectral typing, we needed op- 
tical spectral types obtained with a comparable tech- 
nique, i.e. comparison of equivalent widths between 
the TTS and photometric standards. To this end we 
used archival, low-dispersion optical spectra obtained at 
the 1.5m telescope of the Whipple Observatory with the 
Fast Spectrograph fo r the Tillinghast Telescope (FAST; 
iFabricant etaHfroM ) on the Loral 512x2688 CCD. Spec- 
tral coverage of ^3600-7500 A with a resolution of ~6 
A was achieved with the standard configuration used for 



FAST COMBO projects: a 300 groove mm l grating and 
a 3 wide slit. The data were reduced at the Harvard- 
Smithsonian Center for Astrophysics using software de- 
veloped specifically for FAST COMBO observations and 
were wavelength-calibrated and combined using standard 
IRAF routines. Spectra for all of our targets were ob- 
tained i n 1995 and 1996 a s part of Program 30 (PI: 
Kenyon; iKenyon et al.l 119981 ) and are publicly available 
in the FAST database^. 

3. ANALYSIS: SEPARATING THE EXCESSES AND T TAURI 
PHOTOSPHERES 

To constrain the physical parameters of the regions in 
which the NIR excess is emitted, we need to extract the 
flux-calibrated spectrum of this excess as a function of 
wavelength. Since the excess fills in or 'veils' absorption 
lines, we can measure it by comparing line strengths of 
the program stars with those of photospheric templates, 
in this case the two WTTS. To obtain the absolute flux of 
the excess, we also need the extinction towards the star, 
Ay, which can be estimated by comparing the slopes of 
the program star and the template. However, the line 
strengths of the program star are affected by both veil- 
ing and the star's spectral type (SpT), while the slopes 
are affected by those parameters and Ay, so we need 
to determine SpTs, veiling, and Ay simultaneously from 
the same data set. 

The degree of continuum veiling can be defined by the 
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where the superscript phot indi cates the intrin si c flux 
of the underlying photosphere. lHartigan et al.l (|1995f ) 
use Eq. Q] to fit a single value of r\ to several lines 
at once, and this has become the most common tech- 
niqu e used to determine veiling in the recent literature 
(e.g. lMuzerolle et al.|[2003t iFischer et al.ll2011[ ). Alterna- 
tively, the veiling can be written in terms of the equiv- 
alent width, W\, assuming that the veiling is constant 
over the line: 

W phot (X) 

^(A) 1 + ^ () 



iBasri fc Batalha (1990) compare the two approaches and 
find that veilings determined by the latter method have 
less scatter, as one can reject lines known to experience 
differential veiling or surface gravity effects. We therefore 
use the equivalent width method to determine the veiling 
for each of our targets. 

Equivalent widths of particular atomic lines and molec- 
ular bands are also used to determine SpTs, e.g. 
iHernandez et al.l ((2004). The effects of continuum veil- 
ing on the SpT determination can be eliminated by 
considering instead the ratio of equivalent widths of 
two lines close enough in wavelen gth that thei r veil- 
ing is approximately constant (IBasri fc Batalha 119901 : 
ILuhman fc Riekel H99l IVacca fc Sandelll I2011D . Then 

9 FAST database: http://tdc-www.harvard.edu/cgi-bin/arc/fsearch 
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Fig. 1. — Spectra of the ten stars in our sample. Left panel: The two WTTS stars (V827 Tau and LkCa 3) and the three low-accreting 
stars (V836 Tau, FN Tau, and GO Tau) from 0.8 to 4.5 fim. The 2.8 to 4.5 fim regions of GO Tau and FN Tau, our faintest stars observed 
with LXD, have been smoothed to a lower resolution to increase their continuum signal-to-noise ratio. Right panel: The five higher accreting 
CTTS stars from 0.8 to 4.5 fim, in order of increasing Br7 strength, i.e. accretion rate. 



their ratio should be equal to that of a photospheric tem- 
plate of the appropriate s pectral type. What co nstitutes 
"close enough" is unclear. iHartigan et al.l ()1989f ) find the 
veiling to be constant over 10 to 15A intervals at optical 
wavelength, where the veiling continuum is increasing in 
strength towards the U V. The slope of r\ vs. A beco mes 
flatter in the infrared (jWhite fc Hillenbrand! 120041 ) . so 
the veiling is likely constant over a larger interval at the 
SpeX wavelengths. 

3.1. Continuum Determination and Equivalent Widths 



The first step in the procedure outlined above was to 
measure the equivalent widths of individual lines. This is 
relatively straightforward in certain spectral regions, e.g. 
J or K band, but in z and H band there are many molec- 
ular absorption features overlapping the atomic lines. 
Consequently, even if one does find a "good" absorption 
feature to measure, i.e. one with mainly one contribut- 
ing absorber and good signal to noise, it is still difficult 
to define continuum regions for extracting the equiva- 
lent width. To overcome this issue, we use a technique 
outlined by iBasri fc Batalhal (|1990f ) for identifying the 
continuum wavelengths. To summarize, we choose as 
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continuum points the wavelengths corresponding to the 
highest fluxes in the distribution of fluxes within some 
wavelength bin, after performing a cr-clip to remove obvi- 
ous emission lines. We fit these points with a polynomial, 
divide out the continuum, fit each absorption line with 
a gaussian, and use the analytic expression for the area 
of the gaussian to compute the equivalent width. Uncer- 
tainties are propagated from the original uncertainties in 
the flux, and we employ an additional criterion that the 
equivalent width of the line must be greater than 0.2 A, 
which is the rms uncertainty in regions of the spectra 
without obvious absorption lines. The lines which we 
ultimately decided to use are shown in the continuum 
subtracted and normalized spectra of one of our WTTS 
and two CTTS in Fig. [2] Identifications for each line are 
given in Table [3J 

3.2. Comparison of Trends in W\ vs SpT Between 
WTTS, Dwarfs, and Giants 

We began our analysis with two assumptions. First, a 
WTTS with no infrared excess should be a representative 
photospheric template for a CTTS of the same effective 
temperature; i.e. the WTTS has no intrinsic continuum 
veiling. Second, we assumed the Basri fc Batalha] (1990) 
result that there are optical lines for which dwarf stan- 
dard stars are acceptable representations of WTTS ef- 
fective temperatures, and that if a star is spectral typed 
with these lines, any deviations from the standard dwarf 
trends between W\ and SpT from 0.8 to 2.5 /im are 
due to the effects of surface gravity or differential veiling 
caused by chromospheric emission, as discussed in fJTJ In 
this way, we can use our WTTS standards to identify 
which infrared lines are surface gravity sensitive in TTS. 
Because our particular WTTS are binaries, we operated 
with the additional caveat that the lines we ultimately 
chose needed to be more sensitive to spectral types in the 
range of the primaries, K7 to M2, than to late M types, 
i.e. the spectral types of their companions. 

To confirm that our WTTS standards had no ex- 
cess and to determine qualitatively which of the CTTS 
were the most veiled, we constructed equivalent width 
ys spectral type p lots for the IRTF Spectral Library 
(jRayner et al.ll2009f ) dwarf and giant standard stars over 
a sample of lines, as described in Appendix lAl Since some 
'lines' were obviously spectrally unresolved line blends, 
not all of the features produced a clear trend as a func- 
tion of spectral type. However, there were several lines 
for which the equivalent width depended strongly on the 
spectral type over the range predicted for our targets, ~ 
K5 to M5. 

We then derived optical spectral types for our TTS 
sample from the ancillary optical spectra (described in 
32~4l) with the SPTCLASS tooEJ an IRAF/IDL code 
based on the methods described by Hernandez et al.l 
(2004). The code computes spectral types for low mass 
stars (K to M5) by measuring the equivalent widths 
of 16 spectral features that are sensitive to changes in 
the stellar effective temperature. Each spectral index is 
calibrated using spectroscopic standards observed with 
FAST. For most of the program stars, these spectral 
types agree with those from the literature to within 0.5 
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subclasses. We then overplotted the W\ for our sample 
on the dwarf and giant W\ trend plots, assuming the 
TTS optical spectral types from SPTCLASS. 

From this comparison, there are lines over our whole 
wavelength range for which the WTTS lie on the dwarf 
trend curve at locations consistent with their optical 
spectral type, and in these lines it is easy to see the 
degree of veiling in the CTTS relative to the WTTS (see 
Al I doublet at 1.31270 fjm in Appendix El Fig. [TT] 
for example). In addition to our two WTTS, we iden- 
tified two CTTS that have very low accretion rates and 
therefore are not veiled relative to the dwarf standards 
at z, J, H, or K bands at the resolution of our obser- 
vations: FN Tau and V836 Tau. Both have infrared ex- 
cesses from 5 to iOum indicat ive of dusty disks, as seen 
by Spitzer (jFurlan et alJl2006t ). Because these CTTS are 
single within the limits of our selection criteria (see ij2j 
iKraus et al.l 120111 ) and have little to no veiling at our 
spectral resolution and sensitivity, we can use them as a 
check on the position of the WTTS in W x vs. SpT di- 
agrams to ensure that we are not getting contamination 
from the WTTS' companions in our W\ measurements. 
We will refer to these two stars henceforth as 'weakly 
veiled CTTS' and assume that their SPTCLASS optical 
spectral types are accurate, as we did with the WTTS. 

Using the two WTTS and the two weakly veiled CTTS 
we define rudimentary W\ vs. SpT trends for each line. 
It is obvious from the trends that the five other TTS are 
veiled, and that for most of the lines, this veiling is de- 
generate with the spectral classification. That is to say, a 
K7 star with moderate veiling could have the same equiv- 
alent width as either a K3 or M3 star with little veiling, 
depending on the shape of the trend. In addition, almost 
all of the deep lines (> lA) are surface gravity sensitive, 
as indicated by the discrepancies between the dwarf and 
giant standard trends. This effect is particularly notice- 
able for lines in the H band (e.g. Mg I 1.57721 /im line 
in Appendix [A] Fig. ITT]) . In general, the TTS trends fall 
between the dwarf and giant trends, although the degree 
to which this is true varies on a line-by-line basis. 

There are some lines, e.g. Mg I at 1.1833 /im, that show 
shallower features than can be accounted for by surface 
gravity effects (Appendix lAl Fig. [TTj) . For these lines, 
even the WTTS appear to be veiled. The effect is that 
for these lines, even the WTTS appear to have a later 
infrared spectral type than their optical spectral types. 
This may be an effect of differential veiling. Furthermore, 
for some of the shorter wavelength lines that peak in W\ 
at mid-M types, e.g. the 1.14 /im Na doublet (Appendix 
El Fig. E}, the later three WTTS /weakly- veiled CTTS 
stars lie on the dwarf trend, but our ~K7 WTTS, V827 
Tau, and weakly-veiled CTTS, V836 Tau, have W x that 
are much greater than the K7 point in the dwarf trend, 
making them appear to be later, ~M1 to M2. In the case 
of V827 Tau, the binary companion may be affecting this 
line, but for V836 Tau this may be caused by starspots, 
as it seems to affect some stars more than others of the 
same spectral type. 

To fit the TTS trends, we made both a linear fit to the 
WTTS and weakly-veiled TTS and a reduced \ 2 based 
interpolation between the dwarf and giant trends for each 
line. These fits for a set of typical lines are displayed in 
Appendix [A] along with the fraction of the curve at- 
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Fig. 2. — Continuum normalized SXD spectra for DR Tau (black), the strongest accreter in our sample, DE Tau (dark grey), a moderate 
accreter, and LkCa 3 (light grey), our M2 WTTS over the following bandpasses: (top) i and z bands, (middle-top) z and J bands, (middle 
bottom) H band, and (bottom) K band. Of these regions, H band typically has the best S/N. The absorption features used in this analysis 
are indicated with arrows, with identifications given in Table [3] We note that practica l ly all of the absorption features in these spectra are 
real, spectrally unresolved blends of lines and refer interested readers to Rayncr ct al. (2009) for more detailed information. 



tributed to the dwarf trend vs. the giant trend. 

3.3. Spectral Types 

Since SPTCLASS does not account for continuum veil- 
ing, we computed infrared spectral types for each of the 
veiled CTTS using ratios of equivalent widths of two 
nearby lines to avoid the effects of continuum veiling in 
the measurement of the spectral types. Because the lines 
for which there was no surface gravity dependence were 
typically not 'close enough' in wavelength to use those 
exclusively for such ratios, we ultimately used a com- 
bination of those lines and lines for which the surface 
gravity effects were well-fit by our TTS trend. The de- 
tails of our line selection are given in Appendix [B] and 



in Fig. I12[ we display the line ratios that we ultimately 
used. The resulting spectral types for our sample are 
given in Table 2J 

Comparing the derived infrared spectral types with 
those from the literature and from SPTCLASS, listed in 
Table M a nd H the SPTCLASS SpTs are in agreement 
with the IKenvon k. Hartmannl (J1995D spectral types, 
within the error bars, for all of our program stars ex- 
cept for FN Tau (M3 vs M5) and DR Tau (K3 vs K7). 
FN Tau is highly variable and known to flare, while DR 
Tau is highly veiled, which would make strong metal lines 
looker weaker and earlier. This is likely why our veiling- 
independent infrared spectral type for DR Tau is later 
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(K7-M0). Aside from those stars, our infrared spectral 
types are also consistent within the uncertainties with the 
optica l SpT from SPTCLASS and iKenvon fc Hartmannl 
(1995!). This suggests that the differences between opti- 
cal and infrared spectral classifications in the literature 
may be the result of intrinsic differences between TTS, 
dwarf, and giant standard stars, such as surface gravity, 
that are stronger at infrared than optical wavelengths. 

3.4. Veiling 

Having established the spectral types of our program 
stars, we now use these spectral types to determine veil- 
ings over the 0.8 to 2.3 /zm range. Taking the surface- 
gravity interpolated TTS trends to represent the intrinsic 
equivalent widths for the TTS in each line, we computed 
the veilings, r\, from Equation (J2J|. Initially we included 
all of the lines for which we could determine a TTS trend. 
However, this resulted in a large dispersion in the veiling 
of nearby lines, particularly at J and H band. To explore 
this, we compared these lines across the range of CTTS 
spectral types. Many of the features which produced ex- 
ceptionally large or small values of r\ at z, J, and H in 
the moderately veiled stars were close to lines which, in 
the more highly veiled stars, went into emission. For the 
more moderately veiled stars, these emission lines were 
sometimes too weak to be picked up by our er-clipping 
threshold, which had the effect of both filling in some 
lines, making them appear more veiled, or skewing the 
continuum fits, which made them appear either more or 
less veiled, depending on if the continuum fit too high or 
too low. To produce a sample of veilings that reflected 
the continuum veiling only, we excluded all features at 
wavelengths near emission lines appearing in CI Tau and 
DR Tau and refined our continuum fits for some of the 
0.8 to 1.0 /xm lines to have a longer baseline. The final 
veilings for each wavelength of the sample are presented 
in Fig. [3] The dispersion is significantly less at z band, 
although in H band there is still some scatter. We list 
these final veilings in Table [3] 

3.5. Extinction, Stellar Parameters, and Excesses 

To estimate the amount of extinction along the line 
of sight towards each of our program stars, we use the 
relationship between the observed target fluxes, the (ex- 
tinction corrected) photospheric template fluxes, and the 
veiling of the photospheric template at a given wave- 
length. The degree of veiling should be the same before 
and after extinction correction, so some minor arithmetic 
and the assumption of an extinction curve, A \ /Ay leads 
to: 



Ax 



1 



Ax/A, 



-2.5 log (1 + rx 



FphotW 
Ft,obs{^) 



(3) 



This relationship can be used to compute Ay from 
a linear fit with A\j Ay as the abscissa, 2.51og((l + 
r\ ) F r ,h. n t(\)/Ft._ n hs(\)) as the o r dinate, and Ay as the 
slope (Gullbr ing et alj Il998at ISeperuelo Duarte et al.l 
2008t iFischer et al.ll2011l ). 

To test how much our Ay determination would be 
affected by choosing either one of our binary WTTS 
or a dwarf standard as the photospheric template, we 
construct such a plot for our M2 WTTS, LkCa 3. As 
seen in the top panel of Fig. |U there is a clear change 



in slope (and therefore the implied Ay) between the 
0.6> A x I Ay >0.3 region and the 0.3> Ax/Ay > 0.15 re- 
gion. In the middle panel of Fig. [4) we demonstrate how 
this change in slope affects the determination of Ay . For 
LkCa 3, fitting a line to only the data with Ax/ Ay >0.3 
(the z and y bands) produces Ay=0A5, which is what is 
found for optical determinations for this star. However, 
if one fits a line over the whole range of Ax/Ay, the re- 
sulting slope is steeper than that of the z-y band Ay by 
almost a factor of 2, meaning that we have overestimated 
Ay. Extinction correcting the LkCa 3 spectrum to 0.45 
mag and the M2V standard to it reveals the source of 
the discrepant Ay. there is an intrinsic difference in the 
shape of LkCa3 and M2V standard at H and K bands 
(bottom panel, Fig. [4j. This difference is greater than 
what can be attributed to a change in slope induced by 
the LkCa 3 companion and disappears by 2.5/mi. In con- 
trast, the shape of the excess above the dwarf standard 
photosphere is consistent with the lower surface gravity 
seen in TTS; the WTTS lies between the dwarf and gi- 
ant M2 standards for all wavelengths between 1.4 and 
2.5 pan. 

Given that the dwarf standard continuum shape 
matches the WTTS outside of TiO bands from 0.8 to 
1.3 /xm and that the veilings we measured should be 
least influenced by the binarity of our WTTS over those 
wavelengths, we fit Ay over that region, only, for our 
other program stars. To obtain values for Ax/ Ay over 
our wavelength range w e fit a s plined curve to the Ry 
— 3.1 extinction law of iMathia (|1990T ) and interpolate 
from that curve to the wavelength of each veiling esti- 
mate. Three-sigma uncertainties in Ay are determined 
from the uncertainty in the slope parameter of the linear 
fitting routine and are between 10 and 65% of the total 
values. Excesses above the photosphere are constructed 
by extinction correcting the program star spectra using 
the derived Ay, scaling each dwarf photospheric tem- 
plate by the average veiling at ~1 to 1.3 /im (where the 
veiling is the smallest), and subtracting the scaled tem- 
plates from the TTS. The Ay determinations are plotted 
in Fig. [5j while the values of Ay, with uncertainties, are 
listed in Table |4j and the excesses for the sample are 
plotted in Figs. [Hand [7] as part of the analysis in tj4.11 

After determining the best photospheric template, the 
position of the photosphere with respect to the observed 
spectrum, and the absolute, extinction-corrected flux of 
the photosphere for each our program stars, we compute 
their luminosities and radii. We use the extinction cor- 
rected, scaled flux of the photospheric template at J to 
obtain the absolute J magnitude. Assuming the spec- 
tral types found earlier, which are consistent with dwarf 
spectral types at optical wavelengths, we use the effec- 
tive tempera tures and colors in Table A 5 with Equa- 
tion (Al) in IKenvon fc Hartmannl ([1995) to derive the 
absolute bolometric magnitude for each star, and from 
there the luminosity and radius. Our results are listed 
in columns 4 and 5 of Tabl e HI Finally , we a ssume the 
PMS evolutionary tracks of ISiess e~ al. (2000) to deter- 
mine the mass of each program star and record those val- 
ues in the same table. It should be noted that although 
evolutionary tracks are often used to determine mass, 
this method is highly uncertain and the results depen- 
dent on the tracks assumed. However, the uncertainty is 
systematic, so we do not propagate it here. 
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Fig. 3. — Veiling, r\, as a function of wavelength for the final sample of lines and the six veiled CTTS. Systematic uncertainties are 
estimated as ~0.2 based on the uncertainties in the TTS trends. We note that the actual uncertainties for the more veiled targets, CI Tau 
and DR Tau, are likely larger (~ 0.5-1.0) and mostly dependent on how the continuum is defined. 



Comparing our derived values of Ay and L» with those 
from the literature (see Table [2]) , we note that within 
our uncertainties, we are consistent with previous esti- 
mates in most cases. Our values of Ay are not systemat- 
ically higher or lower than those previously found. Com- 
paring the two targ ets which we have in common with 
iFischer et ail (J2011I ). BP Tau and DR Tau, our Ay for 
BP Tau is lower than their value (0.6 vs 1.75 mag), while 
our Ay for DR Tau is larger (2.0 vs 1.54 mag). For DR 
Tau, their Ay is within the uncertainty of our value, and 
for BP Tau, the difference is clearly the result of fitting 
to different wavelength regions; as seen in their Fig. 7, 
determining Ay over A\j Ay from 0.6 to 0.3 only would 
produce a lower value. 

4. ANALYSIS: PARAMETRIC FITS TO THE EXCESS 

Having derived the relevant stellar parameters, extinc- 
tion along the line of sight, and NIR emission excesses we 
turn our attention to disentangling the contributions to 
the emission excess. From the excess plots constructed 
in 3331 we see that the five most veiled CTTS, DS Tau, 



BP Tau, DE Tau, CI Tau, and DR Tau, have an emis- 
sion excess over the entire 0.8 to ~5 fiva. range, while our 
three weakly veiled CTTS, V836 Tau, FN Tau, and GO 
Tau have no excess at the shorter wavelength end, but 
some small excess from 2 to 3 /im. As discussed in the 
introduction, the NIR excess originates primarily in a 
dust sublimation wall, with potential contributions from 
the longer wavelength end of emission from an accretion 
shock or optically thick gas inside the dust sublimation 
radius. We leave such detailed physical modeling to the 
second paper in this series. Here we opt to derive basic 
quantities such as the characteristic temperatures and 
solid angles for different regions of the excess, which is 
possible due to our absolute flux calibration and the wide 
wavelength coverage of the data. The results of this anal- 
ysis will be used to guide detailed modeling efforts in our 
second paper (Paper II). 



4.1. Characteristic Temperatures and Solid Angles 
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Fig. 4. — (Top panel) Comparison between shape of LkCa 3 and 
an M2V standard. Ay is the slope, assuming no veiling. The 
slope changes around A\/Av=0.3 and again around A\/Ay=0.l8. 
(Middle panel) Plot of the relation given by Equation <£3j) for LkCa 
3 (r\ = 0). Ay determined from the two linear fits to the blue or 
entire wavelengths ranges are listed at top. (Bottom panel) Surface 
gravity dependence of continuum shape: observed LkCa 3 spec- 
trum (red), extinction corrected, and the M2V standard (black) 
and M2 III standard (light grey), scaled to the observed spectrum 
at 1.1 fim. The M2III spectrum is the average of the IRTF library 
M1III and M3III, as the nominal M2III spectrum in the library 
appears closer to an M4 III in shape than an M2 III. Note that al- 
though the continuum shape of the WTTS appears 'later' than the 
dwarf standard of the same spectral type, it lies between the dwarf 
and giant standards and can be better explained by a difference in 
surface gravity. 
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Fig. 5. — Ay determination for the entire sample using Equation 
((3J and dwarf standard stars. Ay is the slope of the linear fit. Fits 
to the whole wavelength range may be skewed by surface gravity 
effects, as demonstrated in Fig. [4] We chose to fit the region least 
affected by surface gravity (0.8 to 1.2 /im; blue line) to obtain our 
final Ays, which are listed in Table(4] 

Because we expect to see the longer wavelength tail 
of blackbody-like emission from an accretion shock, plus 
blackbody-like emission from a dust sublimation wall, 
and a possible third emission component from optically 
thick gas, we decided to fit the excess with a simple 
parametric model consisting of three black bodies at hot, 
warm, and cool temperatures. We were able to achieve a 
set of acceptable fits to the excesses by allowing all three 
temperatures to vary within specified ranges in a grid 
and optimizing the solid angle of each blackbody to best 
fit the combined sum of the blackbodies to the excess, 
identifying the fit that has the lowest value for % 2 per 
degree of freedom. 

The cool temperature component was allowed to vary 
from 500 to 1500 K, with the intention of representing 
a typical dust sublimation temperature. The warm tem- 
perature component ranged from 1500 K to the effec- 
tive temperature of the star, intended to simulate the 
unidentified H band excess (either from more highly re- 
fractory sublimating material or hot, optically thick gas 
just inside the typical dust sublimation region). Both 
the cool and warm components varied in increments of 
100 K. The hot temperature component ranged from 
4000 K up to 8,000 K in increments of 1000 K, and 
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was intended to represent a n accretion shock component 
(|Calvet &: Gullbrinsl I1998D . For temperatures higher 
than ^8,000K, our wavelength range only covers the 
Rayleigh- Jeans tail of the Planck function, so to first 
order the temperature of such a fit is degenerate with 
the solid angle. Therefore when considering the quality 
of the fit, we also required that the flux produced by 
the best-fitting temperature and s olid angle not exceed 
the ob served V band flux listed bv lKenvon fc Hartmannl 
(|1995| ). to reduce this degeneracy . Although these stars 
are variable at V band, the e rror in my is betwe en 0.15 
and at most 0.4 magnitudes (jHerbst et al.l 119941 ). much 
less than the range of error introduced by Ay uncer- 
tainties (0.2 to 1.3 magnitudes). The resulting fits are 
displayed in Figs. [6] and [7] with the temperatures and 
solid angles tabulated in Table O 

We found that none of the TTS were well fit with 
temperatures of 1400 K, typically taken as the dust- 
sublimation temperature. Instead, they required tem- 
peratures consistent with what is expected for the dust 
sublimation temperature of high density, highly refrac- 
tory grains; that is to say in the range of 1600 to 2000 K 
(Hcm lev et al.lll994t iPosch et all 12007( 1. which is consis- 
tent wit h the temperatures of the wall in Herbig AeBe 
systems (|Monnier fc Millan-Gabetll2002f) . An additional 
cool component is required to fit the excesses of DR Tau, 
CI Tau, DS Tau, and FN Tau. The cool component has 
a temperature in the range of 800 K to 1000 K for those 
stars. There are two other groups of stars: V836 Tau 
and FN Tau, which are best fit with a cool component 
and no warm component, and BP Tau, DE Tau, and GO 
Tau, which are best fit by a warm component only; i.e. 
their walls are well-represented by a single-temperature 
blackbody at ^1700 K. In all cases, the solid angles of 
the warm and cool components were much larger than 
the solid angle of the central stars. The hot component 
was the only blackbody whose solid angle was less than 
that of the central star, which is consistent with it arising 
in the accretion shock. The weakly veiled CTTS had so 
little veiling shortward of 1 [xm that the 'excess' in that 
region is completely noise. Consequently, we do not list 
a temperature or solid angle for their hot components. 
The more veiled stars had a range of Thot from 6,000 
K to 8,000 K, c onsistent with the tempera ture range of 
shocks given bv lCalvet fc Gullbrind (|1998l ). 

Since the CTTS DR Tau is also in the sample of 
iFischer et al.l (|2011f l , who fit its excesses with three black 
bodies and conclude that has an H band excess that 
cannot be explained by the shock or wall component, we 
chose to see how well our excesses could be fit using their 
Case A and B temperature sets. In both cases, their hot 
and cool components are fixed at 8000 K and 1400 K, re- 
spectively, which are intended to reflect the temperature 
of a hot shock component and the dust sublimation wall. 
For Case A, the warm component is fixed at 5000 K, in- 
tended to represent a lower energy shock, while in Case 
B it is fixed at 2500 K, intended to represent optically 
thick gas in the inner disk. The temperature sets in both 
Case A and Case B produced poorer fits to the extracted 
excess than our best-fitting values, with x 2 /d.o.f. of 74.2 
and 98.3, respectively, compared with a x 2 /d.o.f. of 10.7 
for our fit. Therefore, while we cannot completely rule 
out the possibility that there is a contribution from op- 



tically thick gas in the inner disk of DR Tau, we do not 
require it to fit our excess. 

4.2. Accretion and Wall Luminosities 

While our spectra do not cover the wavelength region 
over which the excess from an accretion shock would be 
strongest, we can make use of the Br7 emission line at 
2.17/im to estimate the accretion luminosity. By inte- 
grating over the continuum subtracted, flux calibrated 
Br7 line from the excess spectrum, we are able to obtain 
an estimate of the luminosit y in Br7. Then from the 
correlation in Equation (2) of iMuzerolle et all ([1998), 

log(L acc /L G ) = (1.26±0.19)Zo ff (L Br7 /L Q )+(4.43±0.79), 

we use the line luminosity to determine the accretion lu- 
minosity, L acci listed in Table [6l Br7 emission was not 
detected in our two lowest-veiling TTS, V836 Tau and 
FN Tau. GO Tau had a minimal detection. Of the five 
more heavily veiled TTS, DS Tau, BP Tau, and DE Tau 
have accretion luminosities between 10 and 20% of their 
stellar luminosities, while CI Tau and DR Tau are ac- 
creting with 90% and 260% of their stellar luminosities, 
respectively. Comparing the values of L acc found here 
with those from the literature, given in Table [21 we find 
that three of our targets, GO Tau, DS Tau, and BP Tau, 
all have lower L acc by 40 to 60% in our analysis com- 
pared with previous estimates. In contrast, DE Tau, CI 
Tau, and DR Tau all have higher L acc by 25 to 100% in 
this work. Some of the variation is certainly real; for ex- 
ample, GO Tau had no Br7 emission above the noise in 
our spectrum, but clearly has L acc estimates in the liter- 
ature. A portion of the variation in our sample could be 
caused by the differences in our procedure for estimating 
Lsr-y- Since we have flux-calibrated, extinction-corrected 
excesses, we determine Lsr-y directly rather than using 
the equivalent width and an estimate of the continuum 
from photometry from a non-simultaneous observation. 
Ultimately, though, the biggest difference is likely in our 
estimates of Ay via its effect on L„. From L acc , we 
can calculate the mass accretion rates onto the stars, 
M, as M= L acc i?»/(GM,), listed in the fourth column 
of Table U With the exception of DS Tau, all of the 
moderate to strong accrctcrs have accretion rates greater 
than 10 8 Mr,/yr the CTTS average va l ue at 1 to 2 Myr 
(|Gullbring et all Il998al IWhite fc Ghe3 12001 . We can 
also use these independent estimates of L acc to compare 
with the luminosities implied by our parametric fits to 
the excess and to determine the radius of the wall, below 
and in $5] 

As a brief check on the credibility of our fits, we com- 
pare the luminosities of the fitted blackbody components 
with those of the central star and the accretion luminos- 
ity derived from Br7. For each component, we deter- 
mined the luminosity from the temperatures and solid 
angles. These values are reported in Table [6] along with 
the combined luminosity of the system, L to t — L* + L acc , 
and the fraction of the total system luminosity emitted 
by the combined wall components, L wa u = L coo i, L warm , 
or L coo i + L warm where applicable. Since the star and 
accretion shock provide the total energy available to heat 
the physical structure in the circumstellar environment, 
including the wall, it is important that the wall have less 
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Fig. 6. — Parametric fits to the emission excess for V836 Tau, FN Tau, GO Tau, and DS Tau. The observed, extinction-corrected TTS 
spectrum, including both SpeX and Spitzer IRS from Furlan ct al. ( 2006) (black) and WTTS photospheric template (dark grey) are plotted 
as well for reference. The excess above the photosphere (light grey, thick) is rebinned to a lower resolution for display purposes and fit by 
three blackbodies (red) with T^ ot (dashed-dotted), T warm (dashed), and T coo ; (dotted). The combined fit is given by the solid, red line. 
Values for the temperatures and solid angles are given in Table [5] We note that the IRS spectra were not included in the fit, but rather 
plotted for independent comparison. 



luminosity than the total system luminosity; in all of our 
cases, this was true. For all but the highest accreter, DR 
Tau, the wall has between 2 and 10% of the total lumi- 
nosity. For DR Tau it is 17%. We also confirmed that 
the luminosity of the hot component did not exceed L acc 
for any of these stars. In the next section, we discuss the 
implications of these luminosities as well as the charac- 
teristic temperatures and solid angles within the context 
of their physical structure and composition. 



5. DISCUSSION 

5.1. Infrared Spectral Typing 

Based on our comparison of trends between equivalent 
width and spectral type in T Tauri, dwarf, and giant 
stars, we determined the degree to which individual ab- 
sorption lines seen in our TTS were affected by surface 
gravity, finding that in the majority of lines between 1 
and 2.5 /im the T Tauri stars lie at an intermediate (W\, 
SpT) position between the dwarf and giants. This is con- 
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Fig. 7. — Parametric fits to the emission excess for BP Tau, DE Tau, CI Tau, and DR Tau. The meaning of each component is described 
in the caption for Figure [6] 



sistcnt with the results of iLuhman fc Riekd ([1998), who 
found it necessary to interpolate between the dwarf and 
giant trends in order to fit well their TTS at K band. 
In general, for the deeper absorption features the dwarf 
trend has higher equivalent widths than the giants, so if 
one simply uses the dwarf trends to determine spectral 
types for TTS, the TTS will appear to have later spec- 
tral types in the IR than in the optical. Such an effect 
has bee n noted by several authors, e.g. IGullbring et al.l 
(|1998bl) . 

An additional factor contributing to a given star hav- 
ing an apparently later spectral type in the infrared is 
the effect of star spots, for which we see some evidence. 



Many of the lines that peak at late M-types show TTS 
equivalent widths that are larger than predicted given 
their optical spectral type. In some cases, e.g. the 1.25 
/jm K I line, there are no differences between the dwarf 
and giant W\ vs SpT trends over our spectral type range, 
suggesting that the observed effect is not related to dif- 
ferences in surface gravity but instead to the presence of 
cooler spots on the stellar surface. 

We also note that the effects of differential veiling or 
spectrally unresolved contributions from nearby chromo- 
sphcric emission lines considerably complicate the anal- 
ysis of lines shortwards of 1.2 /im. For example, in the 
region surrounding the 1.183 /mi Mg I line (Fig. [5J the 
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Fig. 8. — Probable chromosphcric emission around the 1.183 fim 
feature. The T Tauri stars are V827 Tau (K7.5 WTTS, light grey), 
BP Tau (K7.5 CTTS, dark grey), and DR Tau (K7 CTTS, black). 
Wavelengths corresponding to absorption and emission features are 
indicated by a '+' and labelled with their major absorber. 

highest accreters show emission lines at wavelengths cor- 
responding to C I and Ca I I absorption lines in the solar 
spectrum (See Table 6 of iRavner et~aTll2009D . In DR 
Tau, there is even a suggestion of emission from the Mg 
I line itself, although our spectral resolution is too low 
to confirm this. The result of this, as discussed in t j3.4l is 
that almost all of the TTS have Wi.ig^ consistent with 
M2 to M3 dwarf spectral types, regardless of their optical 
spectral type. Only three of the TTS, GO Tau, FN Tau, 
and LkCa3 have an equivalent width for this feature that 
is consistent with the dwarf trend. This effect, combined 
with the likelihood that star spots contri bute to the Na I 
1.1404 line likely produce the result of iVacca fc SandelJI 
(J2011I) , who use a set of four late-M lines plus the Mg I 
1.1833 line to reclassify TW Hydra from a K7 (optical) 
to M2 (infrared). When we consider their preferred line 
ratio for spectral typing, Na I/Mg I (Fig. [9J, we find 
that all of the CTTS except for GO Tau appear to be 
~M2 by this measure, and the K7 WTTS appear to be 
even later. Whatever the exact causes, this ratio is not 
advisable for determining spectral types. We emphasize 
here how important it is to exercise caution when using 
dwarf standards to interpret TTS line equivalent widths 
for this reason. 

In general, when veiling is taken into account through 
equivalent width ratios, the near-infrared is an ideal 
wavelength range over which to measure spectral types 
for K and M stars. The stellar photospheres are at a 
maximum near z and J bands, and the continuum veil- 
ing due to the accretion shock and/or wall excess is at 
a minimum between z and H bands. We note, how- 
ever, that H band may be superior to J for the highest 
accreting stars, due both to the decrease in differential 
veiling from lines in the accretion shock and because the 
bump in H band due to the lower surface gravity in TTS 
means that late TTS photospheres are almost as bright 
at H as J. Therefore, for more extinguished stars, one 
can achieve better S/N at H . 




1.14040/1.18300, r=-0.90 
J I I I l_ 



K0 



M0 
MK SpT 



M7 



Fig. 9. — Equivalent width ratio between 1.1404 (Na I) and 1.183 
(Mg I) lines, following Vacca & Sandell (2011). Black curve is best- 
fitting ratio trend for the dwarf standards, grey curve is best-fitting 
curve for the giant standards, and red curved is the best-fitting 
curve for TTS (assuming that the TTS are supposed to lie on 
the dwarf trend for Mg I, barring chromosphcric or other effects). 
The red, dashed line represents the linear fit to the TTS trend 
based on where the WTTS actually lie in W\ vs. SpT diagram for 
the 1.183 /im feature. Red crosses are the WTTS, assuming their 
spectral types as determined from optical data, while blue crosses 
are the CTTS with spectral types determined in this analysis. Blue 
dashed lines arc the approximate W\ ratio for TW Hya, as given 
by Vacca & Sandell (2011), after dividing their Na I W\ in half, 
as we only measure the longer wavelength line in the doublet. We 
note that all of the K7 TTS would be assigned spectral types Ml 
to M3 using this diagnostic, based on their ratios of 1.4 to 1.65. 

5.2. Emission Size Scales 

Since we can place only tenuous constraints on the hot 
component properties due to the lack of spectral coverage 
shortwards of 0.8 /im, we simply note that for almost 
all of the accreters the hot component solid angles are 
on the order of 1 to 3% of the stellar solid angle, while 
the luminosity of this component is always less than the 
accretion luminosity derived from Br7 (Tables [5] and [6]). 
Therefore the hot component is not inconsistent with an 
accretion shock, particularly as shock emission is known 
to diverge from a single ^8000 K blackbody spectrum, 
with a significant amount emission in the UV or X-ray 
(|Calvet fc GullbrindU998l: llndebv et al.ll2OO90 . 

In contrast, we can place firm lower limits on the emis- 
sion area of the wall, by making the assumption that the 
wall emission is coming from a grey, optically thick dust. 
In that case, and assuming the wall to be avertical sur- 
face, we can use the total system luminosity and T wa u 
to determine a crude estimate of the dust destruction 
radius: 



Rwall = 



I L^ + L ac 



1/2 



(5) 



Our estimates of the radii of blackbody emission are 
given in units of the stellar radius and in AU in Table 
[7] and are on the order of 0.1 AU or less. We note that 
our radii estimates are slightly smaller than what has 
been found using interferometry for BP Tau, CI Tau, 
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Fig. 10. — Approximate best-fitting scale- factors, £, for the 
height, z = ^H, in terms of the disk pressure scale height, H , as 
a function of inclination angle, i, for the six program stars whose 
excesses fits included a warm blackbody around ~1600 to 1800 K. 
Names of each star are given in plot. Solid circles indicate the 
inclination attributed to the systems, where available. 



and DR Tau (lAkeson et all l2005bt lEisner et all [20071: 
lAkeson et al.ll2005al) . As seen in Table [2j the inner radii 
for these disks is found to be ~0.07 to 0.1 AU, assuming 
a ring model, while our blackbody radii for these disks 
are all 0.06 to 0.08 AU. 

Our wall radii are calculated assuming grey dust, rea- 
sonable for millimeter sized grains. If the grains are 
smaller, their opacity at shorter wavelengths increases, 
and they will reach their sublimation temperature at 
larger radii. One way to test if the grey dust, or black- 
body, assumption is reasonable is to compare the height 
of a wall at the blackbody radius given the solid angle 
found from the parametric fits, £l wa ii, to the expected 
'surface' height of gas at the wall temperature and ra- 
dius, z s = AH. The height, z, of a vertical wall with 
area 2irRz is expressed as z — £,H, where £ is a scaling 
factor and H is the pressure sc ale height of the gas at 
that particular radius, given by (D'Alcss io et alJ li.998): 



H = R 3/2 



kT c 



GM*/j,(T c , p c )mH 



1/2 



(6) 



where T c and p c are the temperature and density of the 
midplane. In our case, R — R wa ii, T c = T wa u, and 
p(T c , p c ) = 2.34 under the assumption that the mid- 
plane is predominately molecular. We calculate H for 
each of the stars and list the results in Table [7J Next, 
we estimate £, which is inclination dependent. Assuming 
the wall is vertical, it can be approximated by a cylin- 
der of radius R wa u centered on t he sta r, as discussed in 
Appendix B of jDullemond et aTl (|2001|) . We use the pre- 
scription given there with a grid of angles between and 
90°to identify the value of £ for which the projected area 
best matches the solid angle of the cool blackbody as a 
function of i (see Fig. ITUl) . For the stars with known 
inclination angles, we overplot the best-fitting £ in Fig. 
[TU1 and record these values in Table 

The best-fitting £ are between 10 and 20. These values 
are unrealistically high, i.e. much greater than that of 



the disk photosphere, typically ~4. This suggests that 
although a composite of blackbody fits is able to repro- 
duce the excess flux, the majority of the emission cannot 
be from truly blackbody grains. If the grains are small 
enough not to have grey opacities, their radii would be 
larger than those given by these blackbody estimates (for 
the same value of the solid angle) and therefore the height 
of the wall required to match the solid angle would ob- 
viously decrease. The requirement that the grains be 
large enough to have approximately blackbody shaped 
emission while small enough not to have grey emission 
should place relatively narrow constraints on the maxi- 
mum grain size in the wall for future physical models. 

5.3. Gas in the Inner Disk vs. a Curved Wall 

Our blackbody radii are less than a factor of two less 
than the radii esti mated from r i ng interferometry models- 
Several authors (lEisner et al.l 120071 : iTannirkulam et all 
120081: lEisner et al.H2009D have suggested that these small 
disk radii indicate that the emission comes from optically 
thick gas inside the wall, based on their assumption of a 
particul ar disk temperature structure and SED/visibility 
fitting. iFischer et al.l (|2011[ ) find an H band excess in 
DR Tau that they fit with component at -2500 to 4000 
K, temperatures which would be consistent with those 
predicted for optically thick gas. While our present data 
cannot be compared with visibilities, in our simple, para- 
metric analysis of the excess we found no evidence for a 
component at ~2500 to 4000 K. In all cases, the best- 
fitting models had a blackbody fitting the bulk of the 
H and K band excess with a temperature <2000 K. 
This temperature is consistent with the most refractory 
grains. 

Given that the temperature of the disk at the dust 
sublimation radius depends on the opacity of the dust in 
the wall, it is not straightforward to say that a particular 
radial location in the disk is 'too hot' for dust without ex- 
ploring the grain properties. Refractory dust (e.g. iron-, 
calcium-, or aluminum-rich silicates) have been found in 
solar-system meteorites. There are suggestions of radial 
gradients of dust composition in these meteorites and 
other solar system bodies (i.e. more iron near Mercury's 
location). A gradient in either grain size or composition 
could produce structure in the wall, with the more re- 
fractory grains could potentially exist closer to the star, 
creating curvature. For at least one target with con- 
tinuum emission from 'inside' the dust-subl i matio n ra- 
dius, follow-up observations by iNajita et al.1 (|2009t ) sug- 
gest that the detected excess is unlikely to be be gaseous, 
but may instead be highly refractory dust. Alternatively, 
curvature in the wall structure could be induced by tak- 
ing into account how the dust sublimation temperature 
changes for different local pressures or how the density 
structure changes wh en d ust settling is included, a s do 
llsella fc Nattal (J2005D and ITannirkulam et al.l (|2007f ). re- 
spectively. 

This option is particularly interesting in light of the 
~800 to 1000 K blackbody component in FN Tau, DS 
Tau, CI Tau, and DR Tau, which required two black- 
bodies to fit the wall emission in our parametric analy- 
sis. If one takes the temperatures of the warm and cool 
blackbodies as characteristic of the same physical struc- 
ture, they may correspond to the inner and outer edges 
of the dust-sublimation front, where the warm compo- 
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nent is in the high density, higher dust sublimation tem- 
perature midplane, and the cool component is in the 
lower densit y, lower dust sublimat ion temperature upper 
layers, a la llsella fc Nattal (|2005| ). On the other hand, 
if one takes into consideration the effects of grain size, 
ou r data appear simil a r to t he excess SEDs computed 
by iTannirkulam et all (|2007f ) for their dust-segregation 
model. 

An additional complication for using SED fitting to 
suggest gaseous emission is the potential for a substantial 
contribution from scattered light from the disk around 1 
to 2 /im, again depending on the grain properties. For 
late K and early M stars, this contribution would have 
a temperature of 3,500 to 4000 K. Full disk models in- 
cluding radiative transfer, grain size distributions, and 
a variety of grain compositions are needed in order to 
determine how large a contribution this is. Such mod- 
els would also enable us to test the radial of the wall in 
a more physically realistic way. We address this in the 
second paper in this series, McClure et al. 2013b (Paper 
II). 

6. CONCLUSIONS 

We have combined several techniques from the litera- 
ture, e.g. veiling independent spectral-typing, determi- 
nation of Ay from veiling, and extracting excess emission 
from veiling, to determine self-consistently basic, model- 
independent properties relating to the star and NIR emis- 
sion excess for a small sample of T Tauri stars in Taurus. 
From this work, our main conclusions are: 



Aside from these conclusions, the analysis here demon- 
strates the importance of obtaining simultaneous, moder- 
ate resolution spectra over a wide span of wavelengths to 
self-consistently determine the properties of both young 
stars and their excesses. It also underscores the impor- 
tance of obtaining high quality spectra of single, weakly- 
accreting T Tauri stars of known optical spectral types at 
infrared wavelengths to use as templates for the classical 
T Tauri stars in one's sample. While in general WTTS 
are similar enough to dwarf standards of the same optical 
spectral type, at infrared wavelengths, there are enough 
differences in surface gravity and chromospheric activ- 
ity that neither dwarfs nor giants are as good as photo- 
spheric templates as a WTTS. As we attempt to unveil 
the excess properties of the inner regions of disks around 
CTTS, it is essential to properly constrain the underlying 
photosphere. 

For now, these results are suggestive but preliminary. 
It is necessary to study a much larger sample of stars 
using physical models of the accretion shock and wall to 
understand the details of the infrared excesses in these 
systems. Doing so is worthwhile, however, to broaden 
our understanding of the energetics and dust evolution 
of the innermost regions of circumstellar disks and how 
they impact the terrestrial planet-forming region. 



The later spectral types and colors often found for 
infrared classifications of T Tauri stars are due to 
differences between the photospheres of TTS and 
dwarf or giant standards, e.g. surface gravity. In- 
terpolation between the dwarf and giant trends to 
fit the TTS allows us to correct for the surface grav- 
ity of the stars to obtain spectral types consistent 
with those in the optical and more consistent veil- 
ings. 

The 0.8 to 5 /xm excess in T Tauri stars can be 
modeled successfully without appealing to emission 
from an inner, gaseous disk. Instead, we find evi- 
dence for emission at temperatures cooler than the 
dust sublimation temperature. 

The additional ~ 800 K blackbody required to fit 
the 3 to 5/xm excesses of three of our program stars 
may be evidence for curvature in their sublimation 
walls or a contribution near 5 \xn\ from the disk 
surface beyond the wall. 

The solid angles of the NIR excess are large, and 
the wall heights required to match them (using a 
blackbody approximation to calculate the wall ra- 
dius) are an order of magnitude too large. There- 
fore the wall must be populated with grains that 
are small enough not to have grey opacities in the 
NIR, e.g. less than tens of microns. 

We explore the latter two conclusions further with 
a physical treatment of the wall and disk in Paper 
II. 
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APPENDIX 
A. TRENDS IN W x VS SPT 



For each line in our, we compare d the equivalent wid ths of the TTS with those of the dwarf and giant standard stars 
in the IRTF spectral type library (|Ravner et al J 12009). The library contains at least one standard star of luminosity 
classes V and III per spectral type between F0 and M9. For each line, we produced a plot of the equivalent width of 
that line as a function of spectral type for both the dwarf and giant standards in the IR TF library. To fit each tre nd, 
we used a non-parametric locally weighted scatterplot smoothing (LOWESS) algorithm (jCleveland fc Devlmlll988fl to 
smooth the data in bins of four spectral subtypes and computed an uncertainty in our fit. Several of the lines for 
which the equivalent width depended strongly on the spectral type are shown in Fig. 1111 

In the panel for each line, we also overplot the location of our two WTTS and two weakly veiled CTTS, including 
error bars on both their W\ and on their spectral type, assuming the spectral types determined from optical data 
using the Hernandez et al. SPTCLASS package. There are three types of trends. In the first, the shape of the WTTS 
as a function of spectral type mirrors that of either the giants or dwarfs. For these lines, the WTTS typically either 
lie directly on the dwarf curve or at values between those of the dwarfs and giants, consistent with the intermediate 
surface gravity of TTS. These lines are typically metal lines with peak equivalent widths in the mid-G to K range. For 
these lines we interpolate between the dwarf and giant trend as x x W\ t dwarf + (1 — x ) x W\. giant = W\ t TTS- Examples 
of these TTS trends are plotted in the first three rows of Fig. Illl and we report 100 x x in the lower right-hand corner 
of each panel of that figure. 

However, for the second type of trend, seen in metal lines that peak at mid-M, e.g. 1.1404 (Na I) and 1.25250 (K 
I), our K7 WTTS and K6.5 weakly- veiled CTTS have equivalent widths consistent with M2 to M5 spectral types. In 
particular, these lines have little or no surface gravity dependence over our spectral type range, and our other WTTS 
and weakly-veiled CTTS do lie on the dwarf curve at the location corresponding to their optical spectral type (Fig. 
[TTl bottom- left two panels). This situation is suggestive of star spots, which have spectral signatures consistent with 
cooler effective temperatures. 

The final trend is when all of the TTS lie below both the dwarf and giant trends, despite strong detections (e.g. Fig. 
[Til bottom-right panel). In this case, there are emission lines near that location that appear in our highest accreters. 
These lines corresponds to transitions seen in the solar spectrum. See §5.11 for more details. 



B. SPECTRAL TYPE DETERMINATION FROM W x RATIOS 

For spectral type determination, we took ratios between all pairs of lines within 0.1 fini of each other, computed the 
WTTS trend of the ratio, and then compared how well the WTTS ratio for a particular line pair correlated linearly 
with SpT over the K5 to M5 range. We then took the ratios with a Pearson correlation coefficient r > 0.8, i.e. strong 
correlations, and we then inspected them by eye. We selected only those ratios for which both the dwarf and giant 
ratio trends showed little scatter and in which the WTTS trend shared the same shape as one or both of the standard 
trends. In the end, we identified six independent pairs of lines that satisfy these criteria, which can be see in Fig. 1121 

We note, however, that even under these circumstances, some of the program stars do not fall on the WTTS trendlinc. 
In particular, DR Tau, our strongest accreter, has emission lines near many of our feature pairs, so we are either unable 
to obtain good fits for some features or for others we may have chromospheric emission in part of the feature. The 
lines and targets for which this occurred are noted in Table [3] 
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Fig. 11. — Top three rows: Trends in W\ of surface gravity sensitive lines as a function of spectral type for our WTTS and weakly 
veiled CTTS (positions indicated by red error bars) and the IRTF spectral library dwarf (black filled circles) and giant (grey filled squares) 
standards. Fits to dwarf and giant trends are given by solid lines, with uncertainties indicated by the hatched fill in the appropriate color. 
The best-fitting interpolation between the dwarf and giant trends for the TTS is also plotted (solid red line). Fraction of the dwarf trend 
contributed to the interpolation is given in the bottom right corner (see text). Bottom row: Trends in W\ of potentially star-spot (1 st and 
2 nd panels) or chromospheric (last panel) lines as a function of spectral type for our WTTS and weakly veiled CTTS. 
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TABLE 1 
Target List 



Name 



Obs. Type 



K s a 
(mag) 



Kmko ' 
(mag) 



Standards 

FS 117 
HD 27761 



photometric 10.05 ± 0.018 
telluric 7.27 ± 0.020 



Program stars 








V827 Tau c 


WTTS 


8.230 ± 0.02 


8.22 


LkCa3 c 


WTTS 


7.423 ± 0.02 


7.41 


V836 Tau 


CTTS 


8.595 ± 0.02 


8.12 


FN Tau 


CTTS 


8.189 ± 0.02 


8.14 


GO Tau 


CTTS 


9.332 ± 0.02 


9.43 


DS Tau 


CTTS 


8.036 ± 0.03 


8.10 


BP Tau 


CTTS 


7.736 ± 0.02 


7.42 


DE Tau 


CTTS 


7.799 ± 0.02 


7.57 


CI Tau 


CTTS 


7.793 ± 0.02 


7.79 


DR Tau 


CTTS 


6.874 ± 0.02 


6.74 



Note. - 
a Ks magnitudes are taken from the 2MASS survey 
HCutri et alj|2003h . 

Kmko measured in this work. Uncertainty on these 
magnitudes is ~0.01 mag. 

c Both WTTS are binaries. At K, LkCa 3 (AB) has a 
separ ation of (Z/48 and a magnitude difference, 5m, of 
0.22 (I White fc Ghedl200H ). At H, V827 Tau (AB) 
has a separation of 0'.'09 and 5m=0.58 magnitudes 
(|Krauset al.|[20TH ). 



TABLE 2 

System parameters from the literature 



Name 



SpT a 



A v a 
(mag) 



L t b 
(Lq) 



T b 



Rwall C 

(AU) 



V827 Tau 


K7 


0.28 


0.79 






LkCa3 


Ml 


0.42 


1.7 






V836 Tau 


K7 


0.59 


0.51 






FN Tau 


M5 


1.35 


0.7 






GO Tau 


M0 


1.18 


0.37 


0.10 




DS Tau 


K5 


0.31 


0.58 


0.21 




BP Tau 


K7 


0.49 


0.93 


0.18 


n ns-?+ 0029 n ioq+o.041 
U.083_ 0041 , 0.123_ 0055 


DE Tau 


AI2 


0.59 


0.87 


0.07 




CI Tau 


K7 


1.77 


0.81 


0.43 


0.097 ± 0.008 


DRTau 


K7 




0.87 


1.01 


0.070 ± 0.026, 0.11 ± 0.041 



Note. — 
a SpT and A v are from IKenvon fc Hartmannl (119951 ). 
b Luminositie s ar e fr om Kcnyon & Hartmann (1995), 
IHartmann et~aTl ([1998) and IMuzerolle et all (|2003t T For V827 
Tau, LkCa 3, V836 Tau, and FN Tau, the luminosities listed are 

bolometric. 

c Radii are fr o m inte rferom etric modeling by | Akeson et al.l (|2005a). 
lAkeson et all (|2005bD , and lEisner et afl l|2007l ). The first value in 
each entry is the radius of a ring model, while the second value (if 
present) is for a uniform disk model. The model for CI Tau was 
assumed to be face-on. 
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White, R. J., & Hillenbrand, L. A. 2004, ApJ, 616, 998 
Zolensky, M. E., Zega, T. J., Yano, H., et al. 2006, Science, 314, 
1735 
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TABLE 3 

Veilings 



A ((im) 


Species 


GO Tau 


DS Tau 


BP Tau 


DE Tau 


CI Tau 


DRTau 


0.81980 


Nal 





d 


0.5 


0.3 


0.1 


2.7 


0.88070 


Nal 





d 


C 


0.3 


C 


C 


0.97880 


Ti I 





0.1 


0.1 


0.3 


0.5 


2,1 


0.99150 


TiO 





0.3 


0,1 


0.3 


C 


C 


1.03440 


Cal 





0.3 


0.1 


0.4 


0.6 


2.8 


1.20870 


Mg I/Si I 





0.1 


0,1 


0,1 


0.6 


2.6 


1.31270 


All 





0.1 


0.5 


0,1 


1.0 


2.5 


1.50270 


Mg I 





0.3 


0.1 


0,-1 


a 


C 


1.57721 


Mg I / Fe I 





0.5 


0,1 


0,1 


0.8 


2.7 


1.67550 


All 





0.2 


0.5 


0.5 


0.7 


2,1 


1.71130 


Mg I 





0.3 


0,1 


0,1 


0.8 


3.4 


1.99340 


Fe I/Si I/Ca I 





0.1 


0.6 


0.7 


a 


3.2 


2.11650 


All 


0.1 


1.1 


0.8 


0.9 


1.7 


6.7 


2.29350 


CO 


0.9 


0.8 


0.8 


1.1 


1.1 


b 



Note. — The formal uncertainties on the non-zero measurements were typically 
0.1, except for CI Tau (typically 0.2-0.4) and DR Tau (typically 0.5-1). 
a Features are in regions of poor telluric correction. 
b Features are below the continuum criterion of 2%. 

c Features are in the vicinity of multiple emission lines, preventing a good continuum 
fit. 
d This line had poor signal-to-noise for this target. 



TABLE 4 
Stellar parameters 



Name SpT 

(optical) 



SpT A v a T e{f L, R, 

(infrared) (mag) (K) (L Q ) (R Q ) 



M* 6 



LkCa3 


M2.0 ± 0.5 




0.5±0.1 


3580 


2 0+"- 1 


c 


C 


V836 Tau 


K6.5 ± 1.5 




1.4±0.5 


4060 


i.3tr« 


2.28 


0.74 


FN Tau 


M3.0 ± 1.0 




1.0±0.4 


3470 


0.8l U n1 


2.46 


0.35 


GO Tau 


M1.0 ± 0.5 


M0±1.5 


1.0±0.6 


3850 


3 +0 - 2 


1.13 


0.59 


DS Tau 


K6.5 ± 1.0 


K7 ±1.5 


1.4±0.5 


4060 


7 +0A 
u -'-0.3 


1.69 


0.77 


BP Tau 


K6.5 ± 1.0 


M0±1.5 


0.6±0.2 


3850 


1 +0 ' 2 


2.28 


0.56 


DE Tau 


Ml. 5 ± 0.5 


Mlil.O 


0.9±0.5 


3720 


n q+ - 5 


2.29 


0.47 


CI Tau 


K5.5 ± 1.0 


K7±1.0 


1.3±0.8 


4060 


6+ ' 7 
u -°-0.3 


1.57 


0.78 


DR Tau 


K3.0 ± 2.0 


M0±1.5 


2. Oil. 3 


3850 


6+ 15 


1.79 


0.57 



Note. — 

a Uncertainties on Ay are 3<t. 

6 As given by the ISiess et all ([20001 1 PMS 

derived T e // and L, as input. 

c Because both WTTS are binary, we do not 



evolutionary tracks, using our 
report stellar masses for them. 



TABLE 5 
Excess parameters 



Name 


■*■ cool 


^^cool 


-L warm 


^ {-warm 


Thot 


^hot 


X 2 /d.cf. 




(K) 


(«.) 


(K) 


(fM 


(K) 


(».) 




V836 Tau a 






1200 


10.9 


8000 


0.005 


4.42 


FN Tau 


900 


33.8 


1700 


1.4 






10.97 


GO Tau 






1650 


1.6 


8000 


0.001 


3.3 


DS Tau 


900 


20.1 


1700 


6.2 


8000 


0.026 


3.82 


BP Tau 






1700 


7.9 


6000 


0.051 


9.81 


DE Tau 






1600 


10.5 


8000 


0.007 


18.19 


CI Tau 


1000 


30.0 


1800 


12.3 


7000 


0.02 


2.49 


DRTau 


800 


277.5 


1700 


52.1 


8000 


0.165 


10.74 



Note. — 
a V836 Tau was not observed with LXD, so the temperature is an 
upper limit and the solid angle is not well constrained. 
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TABLE 6 

Luminosities 



Name 


Lg rl 


L'acc 


M 


^cool 


L"warm 


Lhot 


T b 

(Lq) 


L wa ii/Ltot 


Lhot/Lacc 




(L Q ) 


(^©) 


(Me/yv) 


1^0 J 


(Lq) 


(^©) 


(Lq) 


(Lq) 


V836 Tau a 










0.03 


0.02 


1.97 


0.02 




FN Tau 








0.03 


0.02 




0.79 


0.06 




GO Tau 










0.01 




0.25 


0.04 




DS Tau 


4.3e-5 


0.08 


5.6e-9 


0.01 


0.03 


0.07 


0.78 


0.05 


0.88 


BP Tau 


5.8e-5 


0.12 


1.6e-8 




0.08 


0.09 


1.15 


0.07 


0.75 


DE Tau 


6.8e-5 


0.15 


2.3e-8 




0.08 


0.03 


1.05 


0.08 


0.20 


CI Tau 


1.8e-4 


0.51 


3.3e-8 


0.02 


0.07 


0.03 


1.12 


0.08 


0.06 


DR Tau 


4.7e-4 


1.71 


1.7c-7 


0.08 


0.31 


0.49 


2.35 


0.17 


0.29 



Note. — 
a V836 Tau was not observed with LXD, so the temperature is an upper limit and the solid 
angle is poorly constrained. 

Ltot = L* + L acc 
c L wa ii = L coo i + L warm or L coo i or L war m, depending on which combination of those two 
components are present. 



TABLE 7 
Wall radii and heights 



Name 


J^tvarm 

(ft.) 


r^warm 

(AU) 


H J ix wa rm. 


t b 


V836 Tau 


11.5 


0.121 


0.028 




FN Tau 


4.2 


0.048 


0.031 


11.8 


GO Tau 


5.4 


0.029 


0.016 


11.8 


DS Tau 


5.1 


0.046 


0.018 




BP Tau 


5.1 


0.057 


0.026 


16.3 


DE Tau 


5.8 


0.062 


0.029 


13.4 


CI Tau 


6.2 


0.050 


0.022 


12.5 


DR Tau 


9.9 


0.082 


0.029 


12.7 



Note. — 
a H is the gas pressure scale height at a given 
radius, as defined in Equation ©. 
6 £ = z/H is calculated only for the stars with 
well constrained disk inclinations and warm tem- 
perature component fits. 



